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Abstract 

The eigenphase shifts and mixing parameters for elastic p-d scattering below 
the deuteron-breakup threshold are calculated for the AV14 two-body potential 
using two different numerical methods. The excellent agreement confirms that 
it is possible to perform accurate numerical studies of p-d scattering as well 
as n-d scattering for low energies. The numerical results can be considered as 
benchmarks for future calculations. 
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I. Introduction 



The theoretical investigation of the three-nucleon system is an important tool to 
test our knowledge of the nuclear interaction. This system is simple enough that 
one can obtain accurate solutions for the 3 H and 3 He bound-state systems. Modern 
realistic nucleon-nucleon two-body potential models that fit the data with a \ 2 P er 
datum very close to unity underbind the trinucleon bound states by 0.5 - 1.0 MeV. A 
possible solution to this discrepancy is to include the three-nucleon interaction in the 
model Hamiltonian. Various models of the three-nucleon force have been adjusted so 
that the correct three-body binding energy is obtained. Since the other bound-state 
observables scale with the binding energy, the bound-state system cannot be used to 
test our understanding of the three-body force. 

Initially it was believed that the nucleon-deuteron scattering problem could be 
used to test our understanding of the three-nucleon force; however, it has been shown 
U that most of the scattering data can be reproduced at a good level with only two- 
body interactions. For low-energy scattering the effects of the three-nucleon force 
are usually small; there are nevertheless some discrepancies (such as the nucleon 
analyzing power A y and the deuteron analyzing power iTn) that could be sensitive 
to these interactions. Since much of the experimental data is for p-d scattering at low 
energies, where the Coulomb interaction cannot be neglected, it is necessary to solve 
the scattering equations for this case as well as for the n-d case. 

Various groups have published benchmark calculations showing that it possible to 
obtain accurate solutions of the n-d Faddeev scattering equations for energies above 
the breakup threshold f||, H and using s-wave potentials. Recently, a benchmark 
calculation for n-d scattering has been published at Ei a b = 3 MeV (just below the 
deuteron breakup threshold) using a realistic NN potential || . Conversely, at present 
there has been no thorough study of the accuracy of the solutions for p-d scattering. 
In this paper we present a detailed comparison of the p-d phase shifts and mixing 
parameters obtained by two different configuration-space calculations. These results 
can serve as benchmarks for the low-energy p-d scattering problem. The methods 
used for these calculations can be extended above the deuteron-breakup threshold, 
and some preliminary calculations have been published for this case [Q. 

II. Results 

The two methods used for the results in this paper have been described in previous 
papers. The Pisa group || [7], |8| uses the Pair-correlated Hyper spherical Harmonic 
(PHH) basis to expand the wave function, and the corresponding S-matrix is obtained 
using the complex form of the Kohn variational principle. The Los Alamos-Iowa 
(LA-Iowa) group ||| [ItJ solves the Faddeev-Noyes equations using a spline expansion 



in configuration space with boundary conditions appropriate to two-body Coulomb 
scattering in the asymptotic region. 



For the two-body potential we choose the Argonne AV14 model [TT[, since this 
was used in a previous benchmark paper for n-d scattering |4[]. We consider the 
same parameters and energies in order to show the effects of the Coulomb interaction 
on the results. There is one small difference from the previous calculations, we use 
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h 2 /M (MeV-fm 2 ) = 41.47 instead of the value of 41.473 used previously. The con- 
ventions used for the phase shifts and mixing parameters are the same as in Ref. [|J. 
For completeness, we show in Table 1 the binding energies and scattering lengths for 
this potential. The Pisa values shown in the table were previously published in Ref. 
||, while the LA- Iowa numbers are the results of recent more accurate calculations. 
Also, we give in Table 2 the n-d results for our value of h 2 /M. Our primary results, 
the p-d phase shifts and mixing parameters, are given in Table 3, where one can see 
that the two calculations differ by less than 1%, and in many cases (such as the large 
phase shifts) the differences are much smaller than that. 

III. Conclusions 

The results from two very different numerical methods are in excellent agreement 
for p-d scattering at energies below the deuteron-breakup threshold. This confirms 
that in this case it is possible in practice to perform theoretical studies with the same 
precision as for the n-d problem. These results also provide a benchmark against 
which researchers in the future can test their techniques. 
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Pisa 


LA-Iowa 


£b( 3 H) 


7.684 


7.684 


£s( 3 He) 


7.033 


7.033 




1.189 


1.191 


And 


6.379 


6.376 


2fl pd 


0.941 


0.945 


4/, 
Opd 


13.773 


13.752 



Table 1: Comparison of the binding energies in MeV and scattering lengths in 
fermis determined by the Kohn variational method (Pisa) and the configuration-space 
Faddeev equations (LA-Iowa) for the AV14 NN potential. 
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Pisa 
-0.999 

-17.696 
1.043 



1 MeV 
LA-Iowa 
-0.999 

-17.702 
1.041 



Pisa 
-2.573 

-27.866 
1.205 



2 MeV 
LA-Iowa 
-2.573 

-27.864 
1.203 



Pisa 
-3.903 

-34.822 
1.253 



3 MeV 
LA-Iowa 
-3.901 

-34.812 
1.251 



Si, 



-4.195 
12.384 
3.742 



-4.195 
12.383 
3.742 



-6.654 
20.558 
5.394 



-6.651 
20.551 
5.393 



-7.532 
25.052 
7.256 



-7.526 
25.026 
7.255 



<*§o 

Sin 

2 Z 

Sin 
2 Z 

e 
V 



-47.148 
0.580 

-1.073 
0.651 
0.545 

-0.114 



-47.143 
0.580 

-1.073 
0.650 
0.546 

-0.114 



-61.329 
1.547 

-2.771 
0.716 
1.010 

-0.247 



-61.334 
1.547 

-2.770 
0.716 
1.010 

-0.247 



-70.491 
2.419 

-4.214 
0.780 
1.438 

-0.387 



-70.532 
2.418 

-4.212 
0.780 
1.437 

-0.388 



S3-, 
2 x 

e 
V 



0.125 
-4.139 
14.285 
-1.307 
-0.184 
-1.108 



0.125 
-4.139 
14.287 
-1.309 
-0.185 
-1.110 



0.502 
-6.490 
22.722 
-1.970 
-0.269 
-2.316 



0.503 
-6.486 
22.721 
-1.972 
-0.269 
-2.318 



0.943 
-7.197 
26.398 
-2.761 
-0.258 
-3.804 



0.943 
-7.188 
26.382 
-2.766 

0.258 
-3.807 



01, 
2 Z 

S3n 
2 Z 

e 
V 



-0.015 
0.575 
-1.141 
-0.288 
-0.287 
-0.870 



-0.015 
0.575 
-1.141 
-0.288 
-0.287 
-0.870 



-0.093 
1.528 
-2.974 
-0.309 
-0.522 
-1.573 



-0.093 
1.529 
-2.973 
-0.308 
-0.521 
-1.572 



-0.211 
2.384 
-4.567 
-0.328 
-0.737 
-2.156 



-0.211 
2.384 
-4.564 
-0.327 
-0.735 
-2.154 



% 

Sin 
2 d 

S3n 
2 J 

e 

v 



13.442 
-0.065 

0.131 
-0.468 

0.414 
-0.131 



13.452 
-0.065 

0.131 
-0.468 

0.415 
-0.131 



22.026 
-0.257 
0.523 
0.493 
0.728 
-0.255 



22.029 
-0.257 
0.524 
0.493 
0.729 
-0.255 



26.347 
-0.476 
0.971 
0.517 
0.985 
-0.361 



26.350 
-0.476 
0.971 
0.512 
0.985 
-0.361 



Table 2: Comparison of the n-d phase shifts and mixing parameters (in degrees) de- 
termined by the Kohn variational method (Pisa) and the configuration-space Faddeev 
equations (LA-Iowa) for the AV14 NN potential. 
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Pisa 
-0.787 

-12.607 
1.189 



1 MeV 
LA-Iowa 
^787" 

-12.616 
1.187 



Pisa 
-2.281 

-23.525 
1.257 



2 MeV 
LA-Iowa 
-2.281 

-23.535 
1.256 



Pisa 
-3.615 

-31.400 
1.259 



3 MeV 
LA-Iowa 
-3.615 

-31.414 
1.260 



Si, 



-3.394 
9.431 
3.128 



-3.393 
9.432 
3.133 



-6.122 
17.896 
4.589 



-6.120 
17.894 
4.593 



-7.405 
22.815 
6.227 



-7.401 
22.799 
6.229 



Sin 

2 Z 

Sin 
2 Z 

e 
V 



-37.335 
0.454 

-0.849 
0.819 
0.533 

-0.095 



-37.323 
0.454 

-0.849 
0.817 
0.532 

-0.093 



-53.436 
1.357 

-2.459 
0.792 
0.976 

-0.215 



-53.436 
1.356 

-2.459 
0.796 
0.974 

-0.213 



-63.673 
2.206 

-3.905 
0.833 
1.385 

-0.343 



-63.705 
2.204 

-3.903 
0.836 
1.384 

-0.340 



S3-, 
2 x 

e 
V 



0.100 
-3.358 
10.957 
-1.077 
-0.191 
-0.998 



0.100 
-3.358 
10.958 
-1.080 
-0.190 
-0.994 



0.447 
-6.007 
20.152 
-1.650 
-0.296 
-2.091 



0.447 
-6.005 
20.150 
-1.653 
-0.295 
-2.092 



0.872 
-7.166 
24.628 
-2.327 
-0.322 
-3.362 



0.873 
-7.158 
24.615 
-2.333 

0.321 
-3.365 



01, 
2 Z 

S3n 
2 Z 

e 
V 



-0.011 
0.451 
-0.905 
-0.378 
-0.318 
-0.997 



-0.011 
0.451 
-0.905 
-0.378 
-0.320 
-1.006 



-0.081 
1.342 
-2.640 
-0.353 
-0.530 
-1.617 



-0.081 
1.342 
-2.639 
-0.355 
-0.529 
-1.616 



-0.194 
2.176 
-4.229 
-0.362 
-0.736 
-2.189 



-0.193 
2.175 
-4.227 
-0.364 
-0.734 
-2.186 



% 

Sin 
2 d 

S3n 
2 J 

e 

v 



10.237 
-0.051 

0.103 
-0.234 

0.409 
-0.136 



10.246 
-0.051 

0.103 
-0.223 

0.408 
-0.135 



19.279 
-0.229 
0.464 
0.259 
0.729 
-0.258 



19.289 
-0.228 
0.464 
0.261 
0.729 
-0.257 



24.216 
-0.443 
0.899 
0.360 
0.992 
-0.366 



24.222 
-0.442 
0.899 
0.361 
0.992 
-0.365 



Table 3: Comparison of the p-d phase shifts and mixing parameters (in degrees) de- 
termined by the Kohn variational method (Pisa) and the configuration-space Faddeev 
equations (LA-Iowa) for the AV14 NN potential. 
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